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The energy-momentum tensor and the tensor continuity equation serve as the conservation laws 
of energy, linear momentum, and angular momentum for a continuous flow. Previously, we de¬ 
rived equations of motion for macroscopic electromagnetic fields in a homogeneous linear dielectric 
medium that is draped with a gradient-index antireflection coating (J. Math Phys. 55, 042901 
(2014) ). These results are consistent with the electromagnetic tensor continuity equation in the 
limit that reflections and the accompanying surface forces are negligible thereby satisfying the con¬ 
dition of an unimpeded flow in a thermodynamically closed system. Here, we take the next step and 
derive equations of motion for the macroscopic fields in the limiting case of a piecewise-homogeneous 
simple linear dielectric medium. The presence of radiation surface forces on the interface between 
two different homogeneous linear materials means that the energy-momentum formalism must be 
modified to treat separate homogeneous media in which the fields are connected by boundary con¬ 
ditions at the interfaces. We demonstrate the explicit separation of the total momentum into a 
field component and a material motion component, we derive the radiation pressure that transfers 
momentum from the field to the material, we derive the electromagnetic continuity equations for a 
piecewise homogeneous dielectric, and we provide a lucid reinterpretation of the Jones and Richards 
experiment. 


I. INTRODUCTION 


The energy-momentum tensor is an innate and com¬ 
pelling aspect of energy and momentum conservation in 
a continuous flow pQ. Recently [Ml, we used global 
conservation principles to construct the total energy- 
momentum tensor for a thermodynamically closed sys¬ 
tem consisting of a quasimonochromatic optical pulse and 
a homogeneous simple linear medium that is draped with 
a gradient-index antireflection coating. Regarding the to¬ 
tal energy-momentum tensor and the tensor continuity 
equation as fundamental, we derived equations of motion 
for the macroscopic fields. The formulation of contin¬ 
uum electrodynamics that was derived in our previous 
work EH3 was limited to homogeneous materials with a 
gradient-index antireflection coating. In this article, we 
develop the theory of continuum electrodynamics for the 
more usual situation of a piecewise-homogeneous linear 
dielectric medium. One of the major differences with the 
Maxwell theory is that the Fresnel relations can no longer 
be derived from the application of Stoke’s theorem to the 
Faraday and Maxwell Ampere Laws. Instead, we derive 
the Fresnel relations from the electromagnetic wave equa¬ 
tion and conservation of energy. We obtain the field and 
material components of the total momentum, derive the 
radiation pressure, and derive the electromagnetic con¬ 
tinuity equations for a piecewise homogeneous medium. 
We provide an interpretation of the Jones and Richards 
[8] measurement of the optical force on a mirror immersed 
in a dielectric fluid. 


II. EQUATIONS OF MOTION 


We take as a given that propagation of the electromag¬ 
netic field is characterized by the wave equation 


V x (V x A) + 


d 2 A 

dt 2 


= 0 


( 2 . 1 ) 


in a limit in which absorption can be neglected. Dis¬ 
persion is treated parametrically for the arbitrarily long 
quasimonochromatic fields that are considered here. The 
electromagnetic wave equation, Eq. (2.1), is a mixed 


second-order differential equation that can be written in 
terms of first-order differential equations. To that end, 
we define the macroscopic magnetic field 


B = V x A 


and a second macroscopic field 

T-r n 9A 

= ~c~dt ' 

A Maxwell Ampere-like law 

>xbJ5 = 0 

c ot 


( 2 . 2 ) 


(2.3) 


(2.4) 


results from the substitution of the definitions of the 


macroscopic fi elds, Eqs. (2.2) and (2.3), into the wave 
equation, Eq. ©• Two more equations of motion are 
derived from the definitions of the fields. Thompson’s 
Law 


V • B = 0 


(2.5) 


is obtained from the divergence of Eq. (2.2 


tor identity that the divergence of the cur 


with the vec- 
of any vector 
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is zero. The curl of Eq. (2.3) 


_ — n <9B Vn „ 

v x n-— = — xn 

c at n 


( 2 . 6 ) 


is our variant of the Faraday Law. Taking the divergence 
of Eq. (2.4) and integrating with respect to time, we ob¬ 
tain the Gauss-like law 


v n =-II. (2.7) 

n 

A constant of integration has been suppressed in the ab¬ 
sence of charges. Note that each field equation is alge¬ 
braically equivalent to its counterpart in the macroscopic 
Maxwell equations 


n 2 <9E 

V x B-— =0 

c dt 

( 2 . 8 ) 

V • B = 0 

(2.9) 

1 r9B 

V x E H-— = 0 

c at 

( 2 . 10 ) 

V • (n 2 E) = 0 

( 2 . 11 ) 


Refs. [2HZj. Here, we consider a quasimonochromatic 
optical pulse that passes from one simple linear medium 
into a second such medium through a planar interface at 
normal incidence. In this limit of piecewise-homogeneous 
medi a, w e need only retain the homogeneous parts of 
Eqs. (2.4)—(2.7) to obtain 


nd n 

V x B + —— =0 
c ot 

( 2 . 12 ) 

V ■ B = 0 

(2.13) 

nd B 

V x n- =0 

c at 

(2.14) 

v ■ n = 0, 

(2.15) 


where the fields are connected by boundary conditions at 
the interface between different homogeneous linear ma¬ 
terials. While the derivation of the equations of motion 
for macroscopic fields in a homogeneous medium from 
Eqs. (2.41-(2.71 is obvious, the usage of Eqs. (2.12)- 
(2.15) for piecewise-homogeneous matter remains to be 
investigated. 


III. THE FRESNEL RELATIONS 


if we define II = -nE, which we have every right to do 
under the auspices of Maxwellian continuum electrody¬ 
namics. However, the different sets of motional equa¬ 
tions for macroscopic fields have different tensorial and 
relativistic properties. This means that Maxwellian con¬ 
tinuum electrodynamics, which is fundamentally a vec¬ 
tor theory, admits improper tensor transformations of 
coordinates for the coupled equatio ns [ 7|. The deriva¬ 
tion of the motional equations Eqs. (2.4)—(2.7) from La¬ 


grangian field theory appears elsewhere [7]. Now, it can 
be argued that Eqs. (2.4)—(2.7) violate relativity because 


the Lorentz factor is 'yd = 1 /-y/1 — n 2 u 2 /c 2 instead of 
7 „ = l/y/l — v 2 /c 2 . However, that argument presup¬ 
poses transformations between a coordinate system in 
the dielectric and a coordinate system in a vacuum lab¬ 
oratory (lab) frame. If we consider, instead, transforma¬ 
tions between two coordinate systems in an arbitrarily 
large region of space in which the speed of light is c/n 
[5j, as we s hould , then 7 d is the correct Lorentz fator and 
Eqs. (2.4) (2.7) are, in fact, compatible with relativity 

m- 

Returning to the equations of motion for the macro¬ 
scopic fields, Eqs. (2.41 (2.7), we can identify two lim¬ 


iting cases of particular interest. The first is when the 
gradient of the index of refraction Vn is sufficiently small 
that reflections and Helmholtz forces can be neglected 
[2]. It is in this limit of an unimpeded continuous flow of 
electromagnetic radiation that the energy and momen¬ 
tum densities can be expressed through continuity equa¬ 
tions and a total energy-momentum tensor for a ther¬ 
modynamically closed system. This case was treated in 


First, we demonstrate how our results apply to the 
most obvious issue relating to piecewise-homogeneous 
matter, the Fresnel relations. Because the equations of 
motion for macroscopic fields in a piecewise-homogeneous 
linear medium have changed, there is a question that 
arises as to how the Fresnel relations survive. Apply¬ 
ing Stokes’s theorem to the macroscopic Maxwell curl 
equations, Eqs. (2.8) and (2.10), it was long ago found 


that the tangential components of the electric field E and 
the magnetic auxiliary field H = B are continuous at a 
planar interface between linear homogeneous dielectrics. 
Similarly, the normal components of the displacement 
field D = n 2 E and the magnetic field B were shown to 
be continuous by the divergence th eorem ap plied to the 
Maxwell divergence equations, Eqs. (2.9) and (2.11). The 


simultaneous continuity of the transverse E and H fields 
and the normal parts of the D and B fields leads to the 


Fresnel relations. Because Eqs. (2.6) and (2.7) are inho¬ 


mogeneous, the transverse and normal components of the 
macroscopic field II are not continuous at the material 
interface, even though the appearance of Eqs. (2.14) and 
(2.15) might seem to suggest otherwise. 


We treat the propagation of a quasimonochromatic op¬ 
tical pulse through a piecewise-homogeneous simple lin¬ 
ear medium. The material is initially stationary in the 
laboratory frame of reference and is rigidly attached to 
a support. We take as a given that propagation of the 
electromagnetic field is characterized by the wave equa¬ 
tion 

V«(VxA )t ^=0. (3.1) 
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FIG. 1: Amplitude of the incident field envelope (arb. units) 


FIG. 2: Refracted field entirely within the linear medium and 
separated from the reflected field. 


The pulse is sufficiently monochromatic that dispersion 
can be neglected in accordance with the characterization 
of a simple linear medium. Further, the simple linear 
material with its support is assumed to be arbitrarily 
massive so that n{ r) can be treated as a time-independent 
function of space in the laboratory frame of reference. 
In order to not overly complicate matters, we assume 
normal incidence and adopt the plane-wave limit for the 
field. The vector potential, in the plane-wave limit, 

A (z,t) = ^ (A(z,t)e- i t Udt - kdZ ) + A*(z,t)e^ Ud *- kdg ^ 


this point because their provenance from the macroscopic 


Maxwell equations, Eqs. (2.8 H2.11), is suspect. Because 


the macroscopic field II is not continuous at a step index 
boundary, we present a derivation of the Fresnel relations 
that is based on the wave equation and conservation of 
ener gy. A pplying Stokes’s theorem to the wave equation, 
Eq. (3.1), we have 


(j (V x A) • d\ = / (V x (V x A)) • hda. 

Jr Js 


(3.2) 


can be written in terms of an envelope function A(z, t) 
and a carrier wave with center frequency u>d- Here, kd is 
the amplitude of the wave vector = ( nujd/c)e z that 
is associated with the center frequency of the field, and 
e~ is a unit vector in the direction of propagation along 
the z axis. The vector potential amplitude A is not 
slowly varying if there is a backward propagating field 
component. Figure 1 shows a one-dimensional represen¬ 
tation of the amplitude of the incident field A 0 (z) = 
f ~ - \ 1/2 

A(z,f 0 ) ■ A* (z, to)) about to enter a simple linear 
medium with ri 2 = 1.40 at normal incidence from the 
vacuum ni = 1. Figure 2 presents a time-domain nu¬ 
merical solution of the wave equation at a later time t\ 

depicted by A\(z) = ^A(z,ti) ■ A*(z,ti)j • The re¬ 
flected field and the refracted field have separated and 
the refracted field is entirely inside the medium. The re¬ 
fracted pulse has not propagated as far as it would have 
propagated in the vacuum due to the reduced speed of 
light c/n in the material. In addition, the spatial extent 
of the refracted pulse in the medium is Wt = niWi/ri 2 in 
terms of the width wt of the incident pulse. As shown in 
Fig. 2, the amplitudes of the reflected and refracted fields 
are different from the amplitude of the incident field. 

We would like a formula for determining how much 
of the incident pulse goes into the reflected field and 
how much is refracted. Although the Fresnel relations 
are the necessary formulas, their use is problematic at 


Consider a thin right rectangular box or “Gaussian pill¬ 
box” that straddles the interface between the two medi¬ 
ums with the large surfaces parallel to the interface. 
Then S is the surface of the pillbox, da is an element 
of area on the surface, and n is an outwardly directed 
unit vector normal to da. There is no contribution to 
the surface integral from the large surfaces for our nor¬ 
mally incident field because V x (V x A) is orthogonal 
to n. The contributions from the smaller surfaces can be 
neglected as the box becomes arbitrarily thin. Then 


j (V x A) 


• dl = 0. 


(3.3) 


We choose the closed contour C in the form of a rectan¬ 
gular Stokesian loop with sides that bisect the two large 
surfaces and two of the small surfaces on opposite sides 
of the pillbox. Here, dl is a directed line element that 
lies on the contour, C. Then C , like S , straddles the ma¬ 
terial interface. For normal incidence in the plane-wave 
limit, the field V x A can be oriented along the long 
sides of the contour C. Performing the contour integra¬ 
tion in Eq. (3.31, the contribution from the short sides of 


the contour are neglected as the loop is made vanishingly 
thin and we obtain 


(V x A)i • Ali + (V x A) 2 • Al 2 = 0 (3.4) 


from the long sides, 1 and 2, of the contour. 
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For linearly polarized radiation, we can write the vec¬ 
tor potential of the incident, reflected, refracted, and 
transmitted waves as 

Aj = e x A ie ~ i{uidt - klz) (3.5) 

A r = e x A r e~ i ^ dt+klz) (3.6) 

At = e x A t e- i{u>dt - k2z) (3.7) 


A t = e x A T e~ i{uidt - klz) . 


(3.8) 


where ki = niUd/c , k 2 = muid/c, and e x is a unit polar¬ 
ization vector. It is understood that we use the real part 
of complex fields and neglect double frequency terms in 
field products. For convenience, the scalar amplitudes 
are taken to be real. Using the fact that the line ele¬ 
ments All and AU in Eq. (3.4) are equal and opposite, 
we obtain a relation 


ni(Ai - A r ) = n 2 A t 


(3.9) 


width of the pulse. The refracted pulse has a width of 
niWi/ri 2 due to the change in the velocity of light be¬ 
tween the two media. Then, evaluating the integrals of 
Eq. ( 3.12| ) results in 


n\A 2 = nlA 2 r + n 2 mA 2 t . (3.13) 

Grouping terms of like refractive index, the previous 
equation 


m ( A 2 - A^j = n 2 A 2 (3-14) 

becomes more suggestive as 

n i ^Aj — A r ^ ^Aj + A r ^ = n 2 A 2 (3.15) 


by factoring the binomial. The second-order equation 
can be w ritten as t wo fi rst-order equations. Substituting 
Eq. (3.9) into Eq. (3.15), we have the unique decomposi¬ 
tion 


n 1 



n 2 A t 


(3.16) 


between the amplitudes of the incident, reflected, and 
refracted fields. In order to derive boundary conditions, 
we need another such relation. 

For a stationary simple linear material, the electromag¬ 
netic energy 



(3.17) 


of Eq. ( 3.13 ). We eliminate A t from Eq. (3.16) using 
Eq. (3.17) to obtain 


U= i ^ I w(^) + (V x A) 2 ] dv (3.10) 


A r n\ — ri2 
Ai ni + n 2 


(3.18) 


is conserved. Here, the volume of integration, which in¬ 
cludes all fields present, has been extended to all-space 
a. The total energy is invariant in time by virtue of be¬ 
ing conserved. The total energy at time to, the incident 
energy U(to) = Ui , is equal to the total energy at a later 
time ti, U[t\) = U r + U t , which is the sum of the re¬ 
flected energy U r and the refracted energy U t when the 
refracted field is entirely within the medium. 

In terms of the incident, reflected, and transmitted en¬ 
ergy, the energy balance U(t 0 ) = U{t\) is 


Ui = U r + U t • 


(3.11) 


Substituting Eqs. (3.5) (3.7) into the formula for the 
ener gy, E q. (3.10), and expressing the energy balance, 
Eq. (3.11), in terms of the amplitudes of the incident, 
reflected, and transmitted vector potential results in 


Subsequently, we eliminate A r to get 

A* _ 2m 
Ai ni + n 2 


(3.19) 


We see that the usual Fresnel relations can be derived 
without invoking Maxwell’s equations. Conservation of 
energy, by itself, is sufficient to deri ve Eq . (3.13). How¬ 
ever, there are several ways that Eq. (3.13) can be decom¬ 
posed into two first-order equations. The application of 
Stoke’s theorem to the wave equation guarante es un ique- 
ness of the decomposition represented by Eqs. (3.16) and 

(13.17) and the Fresnel relations, Eqs. (I3.18I) and (13.19). 


IV. MOMENTUM CONSERVATION IN 
PIECEWISE HOMOGENEOUS MEDIA 


n{A~dv = 


nlA 2 dv 


n^Af dv. 


(3.12) 


In order to facilitate the integration of Eq. (3.12), we 


choose the incident pulse to be rectangular with a nom¬ 
inal width of Wi. The pulse has a finite rise time and a 
finite fall time to reduce ringing, but the short transition 
region can be neglected compared to the arbitrarily large 


The correct form for the momentum of the electromag¬ 
netic field in a dielectric is the subject of the century-old 
Abraham-Minkowski controversy HMS1- The currently 
accepted resolution of the controversy, due to Mpller HU, 
Penfield and Haus [18], Pfeifer et al. [TO] , and others, is 
that the issue is undecidable because neither the Abra¬ 
ham momentum nor the Minkowski momentum is the 
total momentum. If we adopt this viewpoint, then the 
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FIG. 3: Transmitted field has left the medium. 


Abraham momentum and the Minkowski momentum are 
irrelevant. Here we use conservation of total energy and 
conservation of total momentum to show that each com¬ 
ponent of the total linear momentum, reflected, refracted 
or transmitted, and kinematic, has a definite expression 
in terms of the macroscopic fields. 

The energy and momentum of an electromagnetic pulse 
in the vacuum 


U v 



(n; + B2) dv 


(4.1) 


G„ = [ BvXUv dv 
Jo C 


(4.2) 


are well-defined and settled. Here, v denotes a quan¬ 
tity that is based in the vacuum, m = 1. The total 
energy and the total momentum of our system in the 
initial configuration at to as shown in Fig. 1 are con¬ 
sidered to be given quantities. Figure 3 shows the re¬ 
sult of continuing the numerical solution of the wave 
equation, Eq. (3.1), until the pulse has propagated com¬ 


pletely through the medium. The incident, reflected, and 
transmitted fields are in vacuum with well-defined ener¬ 
gies so that we can write an energy balance equation 
u,, = U r T Ut T Ukinematic Here, Ukinematic is the kine¬ 
matic energy of a solid block of dielectric material. Writ¬ 
ing the components of the energy balance equation in 
terms of the corresponding vector potential amplitudes, 
we have 






■U, 


kinematic ■ 


(4.3) 


Conservation of linear momentum cause more prob¬ 
lems than conservation of energy because linear momen¬ 
tum is a directed quantity that changes sign upon reflec¬ 
tion. Surface reflection takes momentum from the field 
and transfers the momentum to the material through 
radiation pressure. Once the field has passed entirely 


through the surface, as in Fig. 2, there are no more sur¬ 
face forces and the block of material moves with con¬ 
stant velocity carrying a momentum G kinematic- Note 
that Gkinematic is not the kinetic momentum described 
by Barnett m- Here, in Fig. 3, all the fields have left 
the material and have well-defined momentums in the 
vacuum. Then, we can write 

G i = G r + G T + Gkinematic (4-4) 


or 


I CX.A^G z dv — / QzA r Q z dv 4“ j CxArpG z dv “I- Gkinematic 

J a J <j J a 

(4.5) 

by conservation of linear momentum. Here, a = w^/(2c 3 ) 
is a use ful combi natio n of coefficients. Substituting 
Eq. (4.3) into Eq. (4.5) we find that the kinematic mo¬ 
mentum of the material is 


G kinematic j 2o:A r G z dv Ukinematic /^ ■ (4-6) 

J (J 

Taking Vkinematic! c to be negligible, we find that the 
kinematic momentum is twice the momentum of the re¬ 
flected field, but in the forward direction as determined 
by the direction of the incident field such that 

f _ B r x n, 

G kinematic = 2G r = / 2 dv . (4*7) 

Jo C 

The sign is a bit awkward, but necessary, because B. r x II,. 
is in the backward (negative) direction while the kine¬ 
matic momentum is in the forward direction. 

Now we can return to the situation of Fig. 2 with the 
refracted pulse entirely within the medium. There is a 
linear momentum that is associated with the propagating 
field in the medium [20] although its composition in terms 
of what portion is field momentum and what portion is 
material momentum remains disputed uni. Whatever 
the composition, the momentum that travels with the 
field through the material, G t , must be equal to the well- 
defined momentum of the field 


G t (ti) = G t (0) = [ 

J (7 


Bt X IX7 


dv 


(4.8) 


that has exited the material through the ant ireflec tion 
coating. Applying conservation of energy, Eq. (3.10), we 
find that II t = Villi'. B f = yTiBr, and 


G t (ti) — j 

J a 


B, x IT 


dv 


(4.9) 


is the total momentum that travels with the field inside 
the medium. When comparing Eqs. (4.8) and (4.9), re¬ 
call that the refracted field is spatially nar ro wer than 
the transmitted field. Substituting Eqs (4.7)—(4.9) into 
Eq. (4.41, we find that the total momentum G tot is the 
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sum of well-defined quantities for the refracted momen¬ 
tum, the reflected momentum, and the kinematic mo¬ 
mentum 


Gfnt. — 


B f x n f 


-dv+ 


B r x n, 


-dv— 


b, x n r 


forces are negligible. Then a stationary medium remains 
stationary, the system is thermodynamically closed, and 
energy and momentum are conserved. The tensor conti¬ 
nuity equation is 


-dv. 


J a C J a c J a ( -' 

« 1-Kl! 

The identification of G tot with is proven by demon¬ 
strating conservation of the total momentum. Substi- 


where 


tuting the definitions of the fields, Eqs. (2.31 and (2.4), 


and the Fresnel relations, Eqs. (3.18) and (3.19), into 


d P . = 


d p T afi = 0, 

(5.1) 

n d d d d \ 

(5.2) 

c dtd dx 1 dy 1 dz) 


Eq. ( |4.10[ ) proves that the total momentum is conserved. is the material four-divergence operator [2U51 f7l [2Tj. 


The total momentum has a definite electromagnetic com¬ 
ponent 


^ _ f B t xU t ^ 

G em — / dv 


b, x n r 


dv 


(4.11) 


that is associated with the propagating field. The mate¬ 
rial component 


t~< f o -® r X J 

G ma tl = — / 2- dv 


(4.12) 


is the momentum of the block of dielectric. In the next 
section, we will relate the kinematic movement of the 
block to the Fresnel surface force. Now we consider the 
general case of a dielectric block with index n 2 in an 
inviscid dielectric fluid of index n\. The momentum of 
the field in the dielectric block is known to be 


G em = [ ^^dv 


by Eq. (4.9) and by prior work pHE Then repeal¬ 
ing the above analysis, the total, electromagnetic, and 
material momentums are still given by the formulas, 
Eqs. (4.10) (4.12), although the transmitted and re¬ 
flected fields are different in accordance with the Fresnel 


relations, Eqs. (3.18) and (3.19). 


V. 


ELECTROMAGNETIC CONTINUITY 
EQUATIONS 


T = 


(n 2 + b 2 )/2 (b x n) x 
(Bxn)i w n 

(B x n) 2 w 2 i 

(B X n) 3 W 3 1 


(b x n ) 2 (b x n ) 3 
W 12 W 13 

w 22 w 23 

w 32 w 33 

(5.3) 

is the total energy-momentum tensor PHE and 


Wij = - 11 , 11 , - BtBj 




B 2 )4 


(5.4) 


is the stress tensor PHE 

The imposition of a step-index interface on the incident 
surface of the solid material is accompanied by a surface 
force F due to Fresnel reflection. For a field of cross- 
sectional area A with square temporal dependence, 


AG 


matl 


= -2 A 


B r x n, cAt 
c n\ 


(5.5) 


(4.13) is found by integration of Eq. (4.7). Here, as in the pre¬ 


ceding section, n\ is the refractive index of the region 
from which the field originates, that is, the index of the 
dielectric fluid (or vacuum m = 1) in which the dielectric 
block is immersed. Then 


77-1 AG matl n * B r X IX r 

F =-:-= -2 A - 


At 


(5.6) 


The force must represent a source or sink of electromag¬ 
netic momentum and it must therefore have the same 
time dependence as the momentum continuity equation. 
For a field in the plane-wave limit that is normally inci¬ 
dent on the block of material, the radiation pressure 


The field imparts a surface force to the material due 
to the change of sign of the electromagnetic momentum 
upon reflection. By Newton’s third law, the material ac¬ 
celerates, increasing in momentum. By Newton’s second 
law, the material imposes an equal force on the electro¬ 
magnetic field and momentum is extracted from the field. 
Clearly, the subsystems are open systems as momentum 
is removed from the field and transferred to the material 
by the surface force, but the total system is thermody¬ 
namically closed and the total linear momentum, as well 
as the total energy, is conserved. 

Consider a quasimonochromatic pulse incident on an 
arbitrarily large homogeneous medium. The medium 
is draped with a gradient-index antireflection coating 
and the index changes sufficiently slowly that Helmholtz 


F 

A 


1 Jl\ cAG matl 
A c At 


= —2B r x n r 


(5.7) 


acts on the incident surface at z = 0. Then the radiation 
pressure can be represented in terms of a force density as 


f = (—2B, r x n r )(5(z) • (5.8) 


There is no source or sink of electromagnetic energy so 
we can write a four-force density 

f a = (0, (—2B r x II r )5(*)) . (5.9) 


Then the tensor continuity equation for the unimpeded 
flow of the electromagnetic field, Eq. (5.1), becomes 


dpT als = f a 


(5.10) 
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for a piecewise homogeneous medium. Because the force 
density is a sink of the electromagnetic momentum den¬ 
sity, there is an equal and opposite force that acts as a 
source for the kinematic momentum of the material. Us¬ 
ing the results of Ref. [7] we can derive Newton’s second 
law 


F = M 


dv 

d(t/n i) 


for a material body of mass M immersed in a dielectric 
fluid of index rq. Then the momentum conservation law 
for the solid block of material is 

F = M—^= [ 2B r x U r S{z)dv . (5.11) 

c dt J a 

Note that the tensor continuity equation for a flow of 
non-interacting material particles that is based on a dust 
tensor that is used in Refs. m and [22], for example, 
does not apply here because we have posited a solid di¬ 
electric. The components of the tensor continuity equa¬ 
tion, Eq. (5.101, are the energy continuity equation, 


- — (- (n 2 3 

cdt\ 2 1 


Er 


B x n 


= 0 , 


(5.12) 


and the momentum continuity equation, 
n d 


dt 


;b x n) + v • w = (-2B,. x n r )<5(2). ( 5 . 13 ) 


The momentum continuity equation, Eq. (5.13), explic¬ 


itly displays the relation between the change in momen¬ 
tum on the left-hand side and the effect of the force, 
acting as a momentum sink, on the right-hand side. The 
momentum continuity equation is not an exact compo¬ 
sition of Eqs. (2.12 H 2.151 because boundary conditions 


impose additional constraints. 

One of the enduring questions of the Abraham- 
Minkowski controversy is why the Minkowski momentum 
is so often measured experimentally while the Abraham 
form of momentum seems to be so favored in theoretical 
work. We now have the tools to answer that question. 
The Minkowski momentum is not measured directly, but 
inferred from a measured index dependence of the optical 
force on a mirror placed in a dielectric fluid MU. 
Because the field is completely reflected at the mirror, 
the force on the mirror is 

F = --f-(2cG) = -4- f 2B x n S(z)dv. (5.14) 

c dt c dt Jy 


The measured force on the mirror is directly proportional 
to the refractive index n = n\ of the fluid mm- On the 
other hand, if we were to assume F = dG/dt, then we 
can write Eq. (5.14) as 


F = -4- f 2D x B5(z)dv (5.15) 

c dt J v 


using D = —nil. Then one might infer that the momen¬ 
tum of the field in the dielectric fluid is the Minkowski 
momentum. Instead, we see that the electromagnetic mo¬ 
mentum that is obtained from an experiment that mea¬ 
sures the optical force on a mirror depends on the theory 
that is used to interpret the results. However, based on 
the changes to continuum electrodynamics that are ne¬ 
cessitated by conservation of energy and momentum by 
the propagation of light in a continuous medium, we find 
that Eq. (5.14) is the correct relation between the force 
on the mirror and the momentum of the field in a dielec¬ 
tric. 


VI. CONCLUSION 


The extraordinary persistence of theoretical and ex¬ 
perimental inconsistencies surrounding the Abraham- 
Minkoswski controversy HMS] regarding the energy- 
momentum tensor for light in a linear medium suggested 
the need to re-examine the role of conservation of energy 
and conservation of momentum in classical continuum 
electrodynamics. We found that it is necessary to give 
up the classical macroscopic Maxwell equations in order 
to preserve the tensor form of the energy and momen¬ 
tum conservation laws [Ml- Then we must re-examine 
the body of work that has been built upon the historical 
forms of the macroscopic Maxwell equations. In this arti¬ 
cle, we derived equations of motion, boundary conditions, 
continuity equations and radiation forces for the limiting 
case of macroscopic fields B and II propagating through 
a piecewise-homogeneous linear dielectric medium. 
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